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ABSTRACT 

Popular internal shock models, developed to explain production of high en- 
ergy flares in blazar jets, involve collisions between local overdensities of matter 
being ejected by a central engine and moving along the jet with different veloci- 
ties. Prior to such collisions, the matter is relatively cold and therefore does not 
produce intrinsic non-thermal radiation. However, due to Comptonization of ex- 
ternal radiation by cold electrons, the presence of such matter should be apparent 
by prominent precursor soft X-ray flares, visible prior to non-thermal 7-ray flares. 
In this paper we discuss the predicted properties of such precursors and study 
the dependence of their properties (luminosities and light curves) on kinematic 
parameters of relativistic ejecta and on an angle of view. We demonstrate that 
the lack of evidence for luminous soft X-ray precursors can be reconciled with our 
predictions for their properties if acceleration and collimation of a jet takes about 
three distance decades. We briefly discuss the severe constraints on the internal 
shock models that would be imposed by a non-detection of such precursors. 

Subject headings: galaxies: quasars: general — galaxies: jets — radiation mech- 
anisms: non-thermal — gamma rays: theory — X-rays: general 



INTRODUCTION 



Strong jets in radio loud quasars must be launched very deeply in the gravitational 
potential well, and most likely are powered by rotating black hole and/or innermost parts 
of an accretion disk. Their formation is very likely mediated by magnetic fields, but the 
detailed model is not yet established (compare, e.g., scenarios suggested by Contopoulos 
1995; Lovelace et al. 1997; Vlahakis & Konigl 2003). In magnetic models, acceleration and 
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collimation processes are quite smooth and can proceed over several distance decades. It 
is still not known whether this is accomplished on sub-parsec, parsec, or even larger scales. 
Vlahakis & Konigl (2003) claim that some VLBI data on the variability of the apparent 
structure of parsec scale jets are consistent with the acceleration process in the Poynting 
flux dominated flow. Those authors consider also circular polarization (Homan, Attridge, & 
Wardle 2001; EnBlin 2003) and rotation-measure gradients (Asada et al. 2002; Gabuzda & 
Murray 2003) to be properties that strongly support dynamical dominance of magnetic fields 
in parsec scale jets. However, all these features don't have a unique interpretation: acceler- 
ating radio features, which are actually observed very rarely, can represent shocks formed on 
the interface between a jet and a clump of matter entering the jet from outside and being 
accelerated by the relativistic flow. Likewise, the circular polarization does not necessary 
require dynamically dominated helical magnetic fields (Ruszkowski & Begelman 2002). Fi- 
nally, the gradients in the rotation measure can result from a non-uniform distribution in 
the external Faraday screen (Homan 2004). 

Independent means to verify whether quasar jets are Poynting flux or matter dominated 
is by studying radiative properties of quasars with jets oriented close to the line of sight. 
These quasars together with BL Lac objects form the sub-class of active galactic nuclei 
(AGNs) called blazars. In such objects, any radiation produced within the jet is Doppler 
boosted into the direction of the observer and dominates over other radiative components 
of the AGN. As it was pointed out by Begelman & Sikora (1987) and explored by Sikora & 
Madejski (2000), propagation of relativistic jets through very dense radiation fields present 
in the quasar cores should result in a production of prominent spectral "bumps" in the 
soft X-ray band due to Comptonization of external UV radiation by cold electrons in a 
jet. The lack of such features in the observed blazar spectra indicates that the acceleration 
and collimation of a jet takes place at a distance ~ 10 3 gravitational radii from the central 
supermassive black hole. This strongly supports the idea that jets in quasars are, at least 
initially, magnetically dominated. 

The question whether the jets are still magnetically dominated at distances ~ 0.1 — 3 
pc, where most of blazar radiation is produced, remains open. The most popular blazar 
models invoke internal shocks, assumed to be produced by inhomogeneities moving down 
the jet with different velocities (Sikora et al. 1994; Spada et al. 2001). Since efficient, strong 
shocks cannot be formed in magnetically dominated jets, in these models the energy flux 
of relativistic flow is dominated by matter, specifically by protons. Such jets are likely to 
contain also e + e~ pairs, but with the number n e /n p <tens (as can be inferred from the 
shock energetics; see Appendix A.l and A. 4), they contribute to the inertia of a jet only 
negligibly. Prior to their collision, the inhomogeneities are cold, and, as for the case of the 
cold stationary jet, they upscatter external UV photons up to soft X-ray energies, forming 
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X-ray "precursors" which should precede the actual non-thermal flares. Hence, a comparison 
of the observational data against detailed predictions about the light curves of the soft X-ray 
precursors and the non-thermal flares can provide an efficient tool to search for an evidence 
for inhomogeneities, the basic ingredient of the internal shock models (Madejski, Sikora, & 
Kamae 2003; Sikora & Madejski 2002). If such precursors are not detected, it may indicate 
that non-thermal flares are produced by reconnection events rather than by shocks and, 
therefore, that jets at parsec-scale distances are still dominated by magnetic fields. At the 
very least, this can provide the lower limit to the distance where the jets undergo a conversion 
from the magnetic to matter dominated flows and where shocks could be products of the 
reconnection events. 

In this paper, we develop theoretical framework describing the properties of soft X-ray 
precursors. We show the dependence of their luminosities on kinematical parameters of 
ejecta and on the viewing angle, and illustrate how their light curves should compare with 
the time profiles (including time leads/lags) of the non-thermal flares. 



In our studies of X-ray precursors, we approximate the inhomogeneities by shells which 
propagate ballistically, within boundaries of a cone. In order to illustrate in a simple analyt- 
ical manner the basic relations between the properties of precursors and non-thermal flares, 
we assume that the shells have equal masses and proper widths, denoted as Ao- We also 
assume that the conical angle, 9j, is sufficiently small that the kinematics and dynamics of 
ejecta can be described using a one-dimensional model. With the former assumption, the 
efficiency of the energy dissipation in the colliding shells is maximized, and, therefore, the 
magnitude of the bulk-Compton effects will be somewhat underestimated. Regarding the 
second assumption: the one-dimensional approach can directly be applied for models with 
very narrow jets and, in numerical models, it is generalized for a wider jets by summing 
radiative contributions from the 59 <C 9j shell segments. 



Cold electrons, carried by two neighboring relativistic shells through the external diffuse 
radiation field with their respective bulk Lorentz factors T 2 > Ti 3> 1, upscatter the ambient 
photons with energy hu di ff up to energies 



2. ASSUMPTIONS AND EQUATIONS 



2.1. Bulk-Compton luminosity 



hv BC ~ 



TiVihudiff i = l,2, 
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where 

T>. = _ (0) 

1 r^i-AcosW y> 

and ^ D f, s is the angle between the line of sight and the jet axis. The scattered photons form 
the so-called bulk Compton radiation. That radiation is beamed into the jet direction giving 
the isotropic luminosity 

, dL BC , [*■•<{ 8N e A dE' !c 

where r is the distance at which the jet is fully formed and accelerated; r sh is the distance 
of the discontinuity surface which separates fluids of the two colliding shells; 

dE ' t dE sc i 4 2 

~!T = — = f* c<JTUd ^ (4) 

is the amount of the scattered radiation energy per unit time, where Udi/f is the energy 
density of the external radiation field; and (dN e ^/dr) b s dr is the number of electrons con- 
tributing to the observed radiation at a given instant from a distance range dr. For a uniform 
radial distribution of electrons within the range r^ear < r < rj°^ t 

dN ei \ N e 

9r Jobs W 

and otherwise is zero. Note that the above is due to the fact that the source is observed as 
being stretched by a Doppler factor, i.e. that r^l] t — r^fea) = XoT>i. 

Using Eqs. (3)- (5), one can find that for = const up to a given distance r m and 
then falling very rapidly beyond r m , the maximum bulk Compton luminosity, L BC , is 

4 

L B c,i ^ -caTUdiffT^VfNeei , (6) 

where 

* = Min[l; ^-] . (7) 



2.2. Number of electrons 

The number of cold electrons carried by a pair of shells is equal to the number of 
relativistic electrons involved in the production of non-thermal flares. Assuming that rela- 
tivists electrons are injected with a single power law distribution, Q = K^~ p 1 the number 
of relativistic electrons injected during the collision time scale, t' coll —Vtfi, is 

praax ttlVK 

N e , inj = tflD / Q d 1 = jL- — x = 2N e (8) 
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where tfi is the time scale of the non-thermal flare and p > 1. 

The normalization factor K of the injection function Q can be found from the formula 
for emissivity which in the 7-ray band is very likely dominated by Comptonization of external 
diffuse radiation fields by ultra-relativistic electrons (see, e.g., review by Sikora & Madejski 
2001). At a given frequency we have 

3 V e 

vLv ~ 8 r^ 77 ^ ^ Ecrriec2 ' ( 9 ) 

where iV 7 = dN e /dy is the differential energy distribution of electrons and \^\ec is the ra te 
of the radiative losses of electrons due to the external Compton process. 

Radiation observed in the MeV-GeV band is produced by electrons which cool on time 
scales much shorter than the dynamical time scale (i.e. in the so-called fast cooling regime). 
For them, N 1 ~ J Qd^ /\^\ec, and inserting this into the Eq. (9) gives 

^iS^ 1 ^- (io) 



where 7^ ~ ^Jv E /v di ff /V and v E is the frequency at which the 7-ray luminosity is measured. 

For p — 2, which corresponds with the typical value of the energy flux spectral index in 
the EGRET band a ~ 1 (Pohl et al. 1997), 

" 3m e c 2 P 6 ' { } 



and using Eq. (8) we obtain 



Vt^ 4t fl (u E L UE ) T> 1 

e 9^ ■ r 2 T) 5 ^1 ■ ' 

^ 1 mm <jiiu e Ls ls 1 mm 

where 7 m j„ < a few, as indicated by power-law X-ray spectra extending down to energies 
< IkeV (Tavecchio et al. 2000) and by circular polarization measurements (Wardle et al. 
1998). 

It is worthwhile to note that the number of electrons in cold shells could, in general, 
be much lower than that implied by the non-thermal flares, since a pair production via ab- 
sorption of highest energy 7-rays in the 77 — > e + e~ process is always possible in principle. 
However, the external radiation Compton models predict the 7-ray spectra of quasars having 
the intrinsic cutoff at energies not much larger than 10 GeV (see, e.g. Sikora et al. 1994). This 
together with characteristic 7-ray spectral slopes a ~ 1 indicates that the number of pairs 
produced in the shock cannot significantly exceed the number of the primary electrons and 
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positrons in the ejecta. It also should be mentioned here that in general population of rela- 
tivistic particles in the shocked plasma can be accompanied by population of (quasi-) thermal 
particles. In such a case, the total number of electrons/positrons would be larger and the 
precursors would be more luminous than presented in §3, leading to stronger constraints on 
the internal shock models discussed in §4. 



2.3. The shell width 

In the internal shock model, the time scale of the non-thermal flares is related to the 
time scale of the shock operation. The latter is equal to the time scale of collision which 
depends on the velocities of the shells and on their widths. Since internal shocks are only 
marginally relativistic, the shell width can be derived by using the formula which ignores 
the width of the shocked plasma layer. In the shell rest frame this width is 

A = A' r' 1)2 = ct fl V\(3[^ 2 ~ ctflV^i , (13) 

where |/^ 2 | = (3' 2 = -ft = (a r - l)/(a r + 1), 1% = V 2 = I* = (a r + 1)/(V^)> and 
ar = r 2 /Ti (see Appendix A. 2). Comparison of the approximate value of Ao given by 
Eq. (13) with its exact value is provided in Appendix A. 3. It can be verified using formula 
(A20) that our approximation leads to an overestimate of Ao by no more than a factor 1.4. 
We use the approximate formula in order to calculate self-consistently the light curves, which, 
in addition to the precursors, include the non-thermal flares. The latter are computed by the 
BLAZAR code (Moderski, Sikora, & Blazejowski 2003) which uses thin shell approximation. 



3. RESULTS 

The EGRET instrument detected many prominent flares in blazars (von Montigny et al. 
1995). In the 30 MeV - 10 GeV band, their apparent luminosities reach values vL v ~ 10 48 
erg s _1 , and their time scales, tfi, are on the order of a few days. Relating tfi to the shock 
life-time, one can find that the distance range of the shock operation is 

^ = rrHbu: s ,3ct,,VT K 78 x 1017 {ik) S) 2 ™' (14) 

Roughly symmetrical profiles of the observed non-thermal flares suggest that they occur 
at distances r s h ~ Ar s h (see Sikora et al. 2001). At such distances, the energy density 
of the external diffuse radiation field is dominated by broad emission lines (BEL). From 
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reverberation mapping of the BEL regions in AGNs and quasars we know that these lines 
are produced at 



Tbel ~ 10 18 ( ^ UV I ) cm , (15) 
1 10 4b ergs -1 ' 



0.7 



(Kaspi et al. 2000). 

Energy density of BEL radiation at r < t B el is 

Lbel n „ in -3 /£bel\ ( L uv \ °' 4 _ 3 
MBEL "^V" 2 - 7 10 {-0T){ W« ergs-! J 6rgCm (16) 

where £bel = -^bel/-^uv- At r > t B el the energy density of BEL radiation drops faster 
than r~ 2 and, therefore, in the first approximation, the production of the bulk Compton 
radiation at such distances can be ignored. Another effect which in the first approximation 
can be neglected arises due to the decreasing number of cold electrons due to the collision 
and shock formation. With the above approximations, one can estimate the maximum bulk 
Compton luminosities using Eq. (6) for r m = tbel- The results are presented in Fig. 1, 
which shows the dependence of the precursor luminosities on the bulk Lorentz factor of the 
shocked plasma, and in Fig. 2, which illustrates the dependence of the luminosities on the 
angle of the jet with respect to our line of sight. In both Figures, we plot separately the 
value of L bc for each of the two shells for various values of «r (the ratio of Lorentz factors 
of the two shells); both shells will subsequently undergo a collision, which will result in the 
non-thermal flare. The values of L B q plotted there are calculated by combining Eqs. (6), 
(7), (12), and (13), which give L BC> i oc Y 2 Y 2 Vf /V h for q = 1 (entire source is contributing 
to the observed radiation) and L BC ^ oc T 2 T 2 Vf/V 6 for < 1. 

It is apparent from Fig. 1 that the luminosity of the precursors produced by the faster of 
the two shells becomes much larger than that due to the slower shell when the bulk Lorentz 
factor decreases. The large difference at small T factors reflects the stronger Doppler boosting 
of radiation from the faster moving sources than from the slower ones. The decrease of this 
difference with increasing T is caused mainly by a large reduction of the beaming factor of 
the faster shell for the observer located at fixed angle outside the Doppler cone. 

Fig. 2 shows that bulk luminosities from faster moving shells have lowest values at 
dobs ~ 1/r. For smaller 9 obS) the shape of L BC ,2(9obs) scales roughly with V\\ for larger 9 obs , 
the increase of bulk Compton luminosity is mainly due to increase of a number of electrons 
required to produced a given luminosity of non-thermal flares. Since for 9 b s < 1/r the 
dependence of V 1 on # ofes is very flat, the shape of L BCt i(0 o bs) is determined by a factor V~ 5 
and, therefore, L B q\i increases with the viewing angle. This increase continues at larger 
angles, but is diminished due to the term Vf, which starts to contribute more strongly to 
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the dependence of the L B c,i on the angle of view for 9 obs > l/IY Since the bulk Compton 
luminosities presented on Figs. 1 and 2 are calculated for 9j = 0, while real sub-parsec jets 
are expected to have 9j ~ 1/r, the very high values of L BC2 at 9 Q b s < 1/T are probably 
rather unrealistic. 

The light-curves of soft X-ray precursors and of non-thermal flares are presented in 
Fig. 3 - for jets with 9j <C r 2 , and in Fig. 4 and 5 - for jets with 9j = 9 obs = 1/T. These 
light curves have been calculated for: u^iff oc r -3 for r < = 3 x 10 16 cm; u<nff = Ubel for 
Td < r < Tbel', and ~ r~ 2 ' 5 for r > Tbel- The first branch mimics the contribution of 
the seed photons for the EC process from the accretion disk (Dermer & Schlickeiser 2002). 
This dominates radiative energy losses of relativistic electrons at r < r&. The latter branch 
represents the extension of BEL emissivities beyond the distance r BEL . The time profiles 
of the non-thermal flares for 9j = 1/T are computed using the code BLAZAR described in 
Moderski, Sikora, & Blazejowski (2003). 

It is apparent from Fig. 3 and from the left panels of Figs. 4 and 5 that precursors 
produced by slower shells are very extended in time and reach their maximum luminosities 
significantly before the emergence of non-thermal flares. In contrast, the brighter precursors, 
presumably produced by the faster shells, overlap significantly with the non-thermal flares, 
and, for 9 b s < 1/r, their peaks only marginally precede the peaks of non-thermal flares, 
(see right panels of Figs. 4 and 5). More significant feature of the brighter precursors is much 
faster decay of their light curves than of non-thermal flares. 

In all Figures except for Fig. 3, the cold shells are assumed to be formed at a distance 
rd, above which energy density of external radiation field as measured in a jet comoving 
frame is dominated by the broad emission line radiation. In Fig. 3, the shells are assumed 
to be formed in the close proximity to the central engine and the light-curves plotted in 
those Figures include both Comptonization of broad emission lines and Comptonization 
of the accretion disk radiation. Comptonization of the disk radiation produces additional 
"bumps" in the light curves of the precursors, with a sharp drop at the time corresponding 
to r d . This sudden drop is a result of the apparent lengths of the shells being stretched 
by relativistic effects. For the parameters from Figure 3 the apparent length of the shell is 
much larger than the size of the region in which radiation from the disk dominates, while 
the observed travel time through the region is very short, since it is relativistically reduced 
by a factor (1 — fa cos 9 b s ). This results in a nearly rectangular shape of that part of the 
precursor lightcurve. Figure 3 clearly shows that the Comptonization of the disk radiation 
provides about three times larger luminosity than Comptonization of broad emission lines, 
as indicated by maximum amplitude of the precursors in both cases. 
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4. DISCUSSION AND CONCLUSIONS 

If the 7-ray flares observed in many blazars are produced via internal shocks, they 
should be accompanied by soft X-ray precursors. The precursors are predicted to result 
from Comptonization of external UV continuum and broad emission line photons by cold 
electrons carried by the relativistic ejecta. The ejecta, moving with different velocities, collide 
with each other and form shocks where relativistic electrons are accelerated. The dynamics 
of such events is very complex, and depends on the unknown structure and kinematics of 
the ejecta. However, some basic properties, such as minimal luminosities and profiles of the 
light curves of precursors vs. those of non-thermal flares should not depend strongly on the 
details of the model. We study them by approximating the ejecta by discrete, intrinsically 
identical shells. We analyze the dependence of precursors on various parameters such as the 
bulk Lorentz factor of the shocked plasma, T, the ratio of the bulk Lorentz factors of shells, 
otv = r2/Ti, and the viewing angle, 6Us- 

Our results show that precursors produced by faster shells are typically 10 - 30 times 
less luminous than 7-ray flares if the shells are formed at ro = rj, and about 10 times less 
luminous if the shells are formed much closer to the central engine and contribution from 
the Comptonization of the disk radiation is included. We note here that our model predicts 
that the low energy tail of the EC component of the non-thermal flare contributes in the soft 
X-ray band about 30 - 100 less than at the peak (see, e.g., Moderski, Sikora, & Blazejowski 
2003). The Synchrotron Self-Compton (SSC) contribution is also predicted to be of the same 
order as EC component. With this, the precursors from the faster shells are predicted to 
dominate in the soft X-ray band over the non-thermal flares, even if the contribution from 
the accretion disk is not included. 

As shown in Figures 4 and 5, the brighter precursors overlap significantly in time with 
non-thermal flares, but because their flux should decay much faster than for the non-thermal 
flares, they should be easily distinguished from the synchrotron and SSC contributions, 
provided they are monitored with the intra-day time resolution. Furthermore, presence of 
precursors can be verified by future polarimetric measurements. As predicted by Begelman 
& Sikora (1987), the bulk Compton radiation should be highly polarized, with the electric 
vector perpendicular to the jet axis. The polarized soft X-ray radiation is expected to be 
produced also by synchrotron and SSC process. However, for magnetic fields dominated by 
the transverse component (due to compression of chaotic magnetic fields in the transverse 
shocks), that radiation should be polarized with the electric vector parallel to the jet axis. 
Hence, the polarization position angle measurements and the variability analysis can provide 
two independent diagnosis of the bulk-Compton radiation. 

Clearly, search for the soft X-ray precursors is a viable avenue to study the nature 



of blazar jets, and in particular, to determine if their kinetic luminosities are dominated 
by energy flux of particles or by Poynting flux. The lack of observational evidence for 
precursors by future observations will indicate that cold inhomogeneities originate in situ, 
at distances just prior to the formation of shocks, rather then resulting from modulation of 
the outflow by a central engine. Their formation may be related to the transition from the 
Poynting flux to matter dominated flow; then, the locally fast moving inhomogeneities can be 
a product of intensive magnetic field reconnection events. Another possibility is that the non- 
thermal flares are produced by relativistic electrons directly accelerated in the reconnection 
sites. However, this scenario can be questioned by the fact that the spectra of injected 
electrons predicted theoretically and numerically are much narrower than deduced from 
spectral analysis of the high energy non-thermal flares (Zenitani & Hoshino 2001; Larrabee, 
Lovelace, & Romanowa 2003). In any case, given the promising future of the MeV/GeV 
7-ray astronomy with the upcoming launch of the GLAST mission, this question will be 
addressed directly by observations. However, it is essential that simultaneous observations 
in the soft X-ray band are conducted as well, although the lack of a sensitive all-sky X-ray 
monitor will limit studies to probably a limited number of pointed observations in the X-ray 
band by satellites such as Chandra, XMM-Newton, and Astro- E2. 
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A. APPENDIX 

A.l. Velocity of the shocked plasma and efficiency of the energy dissipation 

Momentum and energy conservation for collision of two intrinsically identical shells gives 
the velocity of the shocked plasma 

= ^^. (AD 

Pi + P2 

and the efficiency of energy dissipation 
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where T = 1/0. - (3 2 . 
For T 2 and Ti > 1 

P - 1 - -wrr , (A3) 

^ 2r!r 2 v ; 

r^^r^f -v^. (A4) 

and 

.. [(r,/r,)v-i]» 

(r,/r,) + i ' 1 ' 

For a given T and a r = r^r^ we have T 1 ~ r/y/oF and T 2 ~ Ty/ar. 

A. 2. Shell velocities in the frame of the shocked plasma 

From Lorentz transformation of shell velocity from the central engine frame to the 
shocked plasma frame 

« = S" (A6) 



For Tj and T ^> 1 it gives 



p2 _ p2 

P'i - ^TF. , (AT) 



r? + r 2 

1 + r : 
2r, ; r 



r 2 + r 2 

H - • (A8) 



For T ~ v / riT 2 the above formulae give 

I" 

2 v^F ' 



r; = r-i^±l, (A9) 



and 

& = -Pi = ■ (A10) 



A. 3. The shock width 

From the shock theory, the ratio of the downstream density na s to the upstream density 
n us is 



r^ = ^ + i (A11) 
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where T' us = T' l 2 is the bulk Lorentz factor of the upstream flow as measured in the down- 
stream (shocked plasma) frame, and 5/3 < 7 < 4/3 (see Blandford & McKee 1976). From 
particle number continuity 

Ao n<is 



. (A12) 

a n us 

where Ao is the proper width of the unshocked shell, and A is the proper width of the 
shocked plasma layer. In the downstream frame (corresponding to that of the shocked 
plasma/discontinuity surface), time of the collision t' coll equals the time necessary for the 
shock front to travel the distance A 

t'coll = -^Ty ( Al 3) 

where f3' s is the velocity of the shock front. In the shocked plasma frame, the width of the 
unshocked plasma shell is \q/T' us and thus during the collision, the outer boundary of the 
unshocked plasma travels the distance 

A /C - A = t' coll c\(3' us \ , (A14) 

where \(3' us \ = \@' 12 \ * s the velocity of the upstream flow in the reference frame of the shocked 
plasma. 



From above equations 



and 



where 



and 



fa - 1)T' 

US ' 

\ = ct' coll \f3>\=ct fl Vg, (A16) 

A = lV ' us + \ = ct fl Vg^ (A17) 
7 - 1 

9 = ^ ~) )l ff us , (A18) 

us ' 

„ i'l^us 1) \Pus \~^us / A 1 n\ 

90 ~ F' + 1 ' ^ ^ 

us ' 



Comparison of the exact value of the shell width given by Eq. (A17) with its approximate 
value given by Eq. (13) gives 

.(approx) A r> -, 

= . (A20) 

^(exact) T"V 1 1 v 1 
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A.4. The pair content 

Amount of energy needed to be injected into relativistic electrons to produce the ob- 
served non-thermal flares is 

/'ymax 
Qmec 2 ^, (A21) 
, min 

where Q is determined from the observed spectra (see § 2.2). This energy cannot exceed the 
amount of energy dissipated during the collision of two shells, which is 

, (Ei + E 2 )r] diss [\oo\ 

^diss = f , { A22 > 

where for n p m p ^> n e m e 

E 1 + E 2 ~ Nprripc 2 ^! + T 2 ) . (A23) 

Hence, for a given amount of electrons determined from observations of non-thermal flares, 
the condition E' diss > E' ein j gives the minimum number of protons in shells, and, re- 
spectively, the maximum pair content. Parameterizing the fraction of dissipated energy 
tapped by relativistic electrons by r) e , and noting that the number of electrons per shell is 
N e = 0.5N e ^ in j = 0.5tfiV J Qd^i we obtain 

jj. = a (Vn7rT- p'_j_ 

where ^ in j = f Q^d^ / f Qdj is the averaged energy of injected electrons. Since n e = n + +n_ 
and rip + n + = n_ (charge neutrality), the pair content (number of pairs per proton) is 
n + /n p = ((N e /N p ) - l)/2. 
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Fig. 1. — Precursor luminosity as a function of the bulk Lorentz factor T. The angle of view- 
is obs = 1/10, and the three pairs of curves are calculated for ratio of the Lorentz factors of 
the two shells ar = 1.25; 2.5; and 5 (marked beside the curves on the left side of the plot). 
The solid lines are for faster of the two shells, the dotted lines are for the slower of the two 
shells. The models are calculated for: u<nff = const = 2.7 x 10~ 3 erg cm -3 for r < tbel and 
u diff — for r > r BEL ; u E L VE = 10 48 erg s" 1 ; tfi = 3 d; and 7 mira = 1. Jet is assumed to 
be very narrow (9j <^ 1/T 2 ). 
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Fig. 2. — Precursor luminosity as a function of the angle of view # {, s , for T = 10. Other 
parameters are the same as in Figure 1. 
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Fig. 3. — Light curves of precursors and flares for four different values of 9 b s . The solid 
lines are the precursors produced by the faster shells; the dotted lines are the precursors 
produced by the slower shells; and the dashed lines are the non-thermal 7-ray flares as 
produced deeply in the fast cooling regime. The models are computed for u^iff = const at 
r d < r < r BEL , and for the power-law distributions at r < r d and r > r BEL as described in 
the text. Precursors starting at r = r d = 3 x 10 16 cm and r = 3 x 10 15 cm are shown by 
the thick and thin curves, respectively. Production of non-thermal flares by shocks starts at 
a distance r 0ySh = 7.8 x 10 17 cm; the shocks terminate at a distance r 0:Sh + Ar sh (see Eq. 14). 
The kinematical parameters are T = 10 and «r = 3. Other parameters are the same as in 
Figures 1 and 2. 
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Fig. 4. — Light-curves of the two types of precursors and the non-thermal flares for o bs — 
(the upper left panel) and 8 f, s = 0.15 (the lower left panel) for 9j — 1/Y — 0.1. Those due 
to the faster shells and the flares are re-drawn in linear scale on the right panels with their 
peaks normalized to one. The models are computed for: tfi — Id; ro yS h = Ar s h, and tq = r<j. 
Other parameters are the same as in Figure 3. 
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Fig. 5. — Comparison of light-curves computed for t/i — 1 day and 3 days. Models are 
computed for 9 obs = 1/T = 0.1. All other parameters are the same as in Figure 4. 



